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An  analytic  model  is  developed  for  scattering  from  random  inhomogeneities  in  range-dependent 
ocean  waveguides  using  the  Rayleigh-Born  approximation  to  Green’s  theorem.  The  expected 
scattered  intensity  depends  on  statistical  moments  of  fractional  changes  in  compressibility  and 
density,  which  scatter  as  monopoles  and  dipoles,  respectively,  and  the  coherence  volume  of  the 
inhomogeneities.  The  model  is  calibrated  for  ocean  bottom  scattering  using  data  acquired  by 
instantaneous  wide-area  ocean  acoustic  waveguide  remote  sensing  (OAWRS)  and  geophysical 
surveys  of  the  ONR  Geoclutter  Program.  The  scattering  strength  of  the  seafloor  on  the  New  Jersey 
shelf,  a  typical  continental  shelf  environment,  is  found  to  depend  on  wave  number  k,  medium 
coherence  volume  Vc,  and  seabed  depth  penetration  factor  Fp  following  a  10  \ogw(F pV ckf) 
dependence.  A  computationally  efficient  numerical  approach  is  developed  to  rapidly  compute 
bottom  reverberation  over  wide  areas  using  the  parabolic  equation  by  exploiting  correlation  between 
monopole  and  dipole  scattering  terms  and  introducing  seafloor  depth  penetration  factors.  An 
approach  is  also  developed  for  distinguishing  moving  clutter  from  statistically  stationary 
background  reverberation  by  tracking  temporal  and  spatial  fluctuations  in  OAWRS  intensity 
images.  ©  2008  Acoustical  Society  of  America.  [DOI;  10.1121/1.2832509] 
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I.  INTRODUCTION 

In  this  paper,  we  develop  an  analytic  model  with  nu¬ 
merical  implementation  for  reverberation  in  range-dependent 
ocean  environments  and  bistatic  source-receiver  geometries. 
The  model  is  derived  from  first  principles  using  the  acoustic 
wave  equation  for  inhomogeneous  media  and  by  application 
of  Green’s  theorem. 

The  model  takes  into  account  the  full  three-dimensional 
(3D)  scattering  interaction  of  the  acoustic  wavefield  with 
volume  or  surface  inhomogeneities.  For  efficiency,  the  model 
is  implemented  in  terms  of  scattering  from  the  spatially  vary¬ 
ing  resolution  footprint  of  the  sonar,  typically  determined  by 
beamforming  and  temporal  matched  filtering.  Statistical  mo¬ 
ments  of  the  randomly  scattered  field  are  expressed  in  terms 
of  statistical  moments  of  fractional  changes  in  medium  com¬ 
pressibility  and  density,  and  the  waveguide  Green  function. 
Variations  in  compressibility  and  density  lead  to  monopole 
and  dipole  terms,  respectively.  As  the  waveguide  Green  func¬ 
tion  is  employed,  the  model  directly  incorporates  attenuation 
in  the  sea  bottom  and  scattering  of  evanescent  waves  by 
inhomogeneities. 

The  model  is  calibrated  for  ocean  bottom  reverberation 
with  data  acquired  by  instantaneous  ocean  acoustic  wave¬ 
guide  remote  sensing  (OAWRS)1"  during  the  2003  Main 
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Acoustic  Experiment  and  geophysical  surveys  of  the  ONR 
Geoclutter  Program.  The  data  were  collected  in  the  New  Jer- 
sey  Strataform,  south  of  Long  Island,  NY,  a  typical  conti¬ 
nental  shelf  environment  with  minimal  bathymetric  relief. 
The  model  is  integrated  to  bathymetric  databases  to  enable 
charting  of  predicted  reverberation  in  geographic  space.  An 
approach  is  also  developed  for  distinguishing  moving  clutter 
from  statistically  stationary  background  reverberation  by 
tracking  temporal  and  spatial  fluctuations  in  OAWRS  inten¬ 
sity  images. 

Analysis  with  the  theory  and  data  indicates  that  (a)  both 
random  fluctuations  in  seabed  compressibility  and  density 
are  important  contributors  to  the  scattered  intensity  and  con¬ 
sequently  reverberation,  (b)  the  scattering  strength  of  the  sea¬ 
floor  on  the  New  Jersey  shelf,  a  typical  continental  shelf 
environment,  is  found  to  depend  on  wave  number  k,  medium 
coherence  volume  Vc,  and  seabed  depth  penetration  factor  Fp 
following  a  10  log1o(ApVcA:4)  dependence,  where  each  term  is 
a  function  of  frequency. 

A  computationally  efficient  numerical  approach  is  devel¬ 
oped  here  to  rapidly  compute  bottom  reverberation  over 
wide  areas  using  the  parabolic  equation4  by  exploiting  cor¬ 
relation  between  monopole  and  dipole  scattering  terms  and 
introducing  seafloor  depth  penetration  factors.  Without  inno¬ 
vations  of  this  kind,  bottom  reverberation  models  based  on 


0001  -4966/2008/1 23(3)/1 270/1 2/$23.00  ©  2008  Acoustical  Society  of  America 


Downloaded  21  Dec  201 1  to  18.38.0.166.  Redistribution  subject  to  ASA  license  or  copyright;  see  http://asadl.org/journals/doc/ASALIB-home/info/terms.jsp 


volume  scattering  can  be  computationally  intensive  in  range- 
dependent  environments  with  bistatic  source-receiver  geom¬ 
etries  because  dipole  scattering  involves  multidimensional 
spatial  derivatives  of  the  complex  waveguide  Green  function. 
These  typically  must  be  computed  numerically  using  large 
and  dense  3D  matrices. 

A  number  of  half-space  models  for  seafloor  scattering 

have  used  the  Rayleigh-Born  approach  for  local  waterborne 

^—8 

backscattering  investigations.  Some  monostatic  range- 
independent  models  have  used  approximations  similar  to 
Rayleigh-Born.1  Previous  range-dependent  waveguide  re¬ 
verberation  models  using  the  parabolic  equation  have  only 
included  monopole  scattering  with  empirical  models.1”-14 


II.  ANALYTIC  FORMULATION 

Here  we  develop  an  analytic  model  for  the  scattered 
field  from  random  medium  inhomogeneities  by  application 
of  Green’s  theorem.  Assume  the  medium  compressibility  and 
density  vary  randomly  in  space,  following  a  stationary  ran¬ 
dom  process  within  the  local  sonar  resolution  footprint. 

Let  the  origin  of  the  coordinate  system  be  placed  at  the 
air-water  interface  with  the  positive  z-axis  pointing  down¬ 
ward.  Let  the  coordinates  be  defined  by  r0  =  (x0,y0,Zo)  f°r  the 
source,  r  =  (x,y,z)  for  the  receiver,  and  r,  =  (xf,yf,z,)  for  the 
centroid  of  an  inhomogeneity.  Spatial  cylindrical 
and  spherical  systems  ( r,0,<p )  are  defined  by  x 
=  r  sin  0cos  </>,  y  =  r  sin  0  sin  <^>,  z  =  r  cos  0  and  p2=x2+y2. 
The  acoustic  wave  number  k  is  given  by  the  angular  fre¬ 
quency  oj = 2  irf  divided  by  the  sound  speed  c. 


A.  Time  harmonic  scattered  field  moments 

In  the  presence  of  volume  inhomogeneities,  the  time- 
harmonic  acoustic  field  satisfies  an  inhomogeneous 

Helmholtz  equation,13 

V2^t(r„,f)+k2^,(rt,f) 

=  -  k2rK(r,)nr„f)  ~  V  •  [T,(r,)  V  <L,(r„/)],  (1) 

where  T(<.  is  the  fractional  change  in  compressibility. 


K 

given  compressibility  K=lldc2,  and  Td  is  fractional  change 
in  density. 


d(  r,)  -  d 
d(  r,) 


(3) 


where  k  and  d  are  the  mean  compressibility  and  density  in 
the  region,  respectively. 

Then,  by  application  of  Green’s  theorem,  given  a  source 
at  r0  and  a  receiver  at  r,  the  time-harmonic  scattered  field 
<5,s(rJ|r,r0,/)  from  inhomogeneities  within  the  sonar  reso¬ 
lution  footprint  Vs  centered  at  rs  is15 


^(rJr»ro,/)  =  J  J  J  [A:2rK(r,)cLf(r„/)G(r|r(,/) 

+  r/Tt)  V  3>,(r t,f)  •  VG{r\r„f)\dVt.  (4) 

Here,  <I>,(ro/)  is  the  total  field  at  the  location  of  the  inhomo¬ 
geneity,  and  G(r|r,,/)  is  the  medium’s  Green  function.  El¬ 
emental  spatial  variations  in  compressibility  lead  to  mono¬ 
pole  scattering,  while  those  in  density  lead  to  dipole 
scattering.  Density  fluctuations  contribute  to  both  monopole 
and  dipole  scattering  while  sound  speed  fluctuations  contrib¬ 
ute  only  to  monopole  scattering.13  The  dipole  term  depends 
on  spatial  gradients  of  the  total  field  and  those  of  the  Green 
function. 

The  total  field  at  any  location  is  a  sum  of  the  incident 
and  the  locally  scattered  fields, 

,  (5) 

where  ch,(r,|r0,/)  =  (4Tr)2G(rf|r0,/).  The  (47t)2  normaliza¬ 
tion  factor  yields  the  convenient  source  level  of  0  dB  re 
1  /rPa  at  1  m.  For  small  local  perturbations  in  medium  com¬ 
pressibility  and  density,  the  total  field  at  the  inhomogeneity 
can  be  approximated  by  the  incident  field  in  Eqs.  (5)  and  (4). 
This  is  the  first-order  Rayleigh-Born  approximation  to 
Green’s  theorem.  The  randomly  scattered  field  at  the  receiver 
can  then  be  expressed  as 

^(rjr,r0,/) 

=  47 T2  III  [£2r*(r,)  G(r,|  r0,f)G(r\rt,J) 

^ s 

+  Td(rt)  V  G(r,|r0)  •  VG(r\rt,f)]dV„  (6) 

which  is  a  single-scatter  approximation. 

The  mean  scattered  field  from  the  inhomogeneities  is 

<^(rJr,r0  ,/)> 

^iJL  [^2<rK(r,))(G(r,|r0,/)G(r]r„/)) 

^ s 

+  <r>,)><VG(rf|r0,/)  •  VG(r\r„f))\dVt.  (7) 

Random  variables  and  T d,  are  assumed  to  be  independent 
of  possible  random  fluctuations  in  the  waveguide  Green 
function.16 

The  second  moment  of  the  scattered  field  is 

<|<bj(rjr,r0,/)|2> 

=  kJh'OJ  L  r0,/)c&*(rs.|  r,  r0,/)> 

=  ((47 r)4  IfL  [/c2T  K(r,)  G(r,|  r0,/)  G(r  |  r„/) 

+  rfi(r,)  V  G(r,|r0,/)  •  VG(r\r„f)]dV, 

Xl  I  /y,^rkr;)G*(r;|r0,t)G*(r|r;,/) 

+  r,(r;)  V  G*(r;|r0,/)  •  VG*(r\r;,f)]dV't).  (8) 

Equation  (8)  requires  evaluation  at  two  spatial  locations  rf 
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and  r',  which  requires  knowledge  of  the  spatial  correlation 
function  of  volume  inhomogeneities.  A  variety  of  different 
forms  have  been  assumed  for  this  correlation  function  in  the 
case  of  seabed  inhomogenieties.  They  include  Gaussian,8 
exponential,  as  well  as  power  law  ’  ’  correlation  func¬ 
tions. 

To  capture  the  essential  physics  in  a  manner  conducive 
to  analysis,  we  assume  that  inhomogeneities  are  uncorrelated 
when  the  two  points  lie  outside  the  coherence  volume  Vc  and 
are  fully  correlated  within  this  volume.  This  leads,  for  ex¬ 
ample,  to 

(T,(r,)r,(r;)) 

=  Vc(rs,zf)[(r^(r,))-|<ri.(rt))|2]^(r,-r;) 

+  <r*(r,)><ri(r;)> 

=  ^c(rs,2f)Var(rit(r,))^(r,  -  r,')  +  (Tk(rt))(Tk(r't)). 

(9) 

Such  delta  function  correlations  are  often  used  in  theories  of 
wave  propagation  in  random  media.1618  The  coherence  vol¬ 
ume  can  be  a  function  of  horizontal  position  and  depth.  For 
example,  overburden  pressure  in  sediments  leads  to  greater 
compaction  with  depth  which  affects  coherence  volume.19 

Substituting  equations  similar  to  Eq.  (9)  into  Eq.  (8), 
leads  to  the  second  moment  of  the  scattered  field, 

<|<Fs(r,|r,r0,/)|2) 

=  |(cFi(rJ|r,r0,/)>|2  +  Var(<f  s(rs|r,r0,/)) ,  (10) 

proportional  to  the  total  intensity,  which  is  a  sum  of  coherent 
and  incoherent  components.  The  relative  importance  of  each 
of  these  terms  within  the  sonar  resolution  footprint  can  be 
estimated  by  analysis  similar  to  that  in  Appendix  A  of  Ref. 
16. 

The  coherent  term  is  the  square  of  the  mean  scattered 
field,  the  square  of  Eq.  (7).  The  incoherent  term  is  the  vari¬ 
ance  of  the  scattered  field, 

Var(<t>s(rj|r,ro,/)) 

=  (4t t)4  J  J  J  Vc(rs,zt) 

V s 

X[k4Var(ri(rf))(|G(rf|r0,/)|2|G(r|r?,/)|2) 

+  Var(rrf(rf)K|VG(r,|r0,/)  *  VG(r|r„/)|2) 

+  k2Cov(r,,rd)(2in{G(r,|r0,/)G(rlr„/) 

X  V  G*(r,|r0,/)  •  VG*(r|r„/)})]<7Vf.  (11) 

Measured  reverberation  from  the  sea  bottom  is  often  domi¬ 
nated  by  the  incoherent  term  of  Eq.  (10).  Exceptions  occur, 
for  example,  when  many  scatterers  with  large  impedence 
contrast,  such  as  layers  of  gravel  or  rock  on  a  sandy  seafloor, 
fall  within  the  sonar  resolution  footprint. 

The  approach  presented  may  also  be  applied  to  model 
scattering  from  irregular  interfaces  between  two  media  with 
differing  sound  speeds  and  densities.  The  volume  integrals  in 
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Eqs.  (7)  and  (11)  for  this  case  would  be  replaced  with  surface 
integrals  over  the  interface  and  the  coherence  volume  Vc 
replaced  by  the  coherence  area  Ac.20 

B.  Broadband  scattered  field  moments 

For  a  potentially  broadband  source  function,  analytic  ex¬ 
pressions  are  provided  for  the  statistical  moments  of  the  scat¬ 
tered  field  from  inhomogeneities  within  the  sonar  resolution 
footprint.  Let  q{t)  be  the  source  wave  form  with  Fourier 
transform  Qif) .  The  expected  time-dependent  scattered  field 
^(rs|r,r0,t)  is 

r  oo 

<%(r,|r,r0,f)>=  G(/)(<Fi(rI|r,r0,/))e^2^/,  (12) 

J  —CO 

by  inverse  Fourier  transform  of  Eq.  (7),  weighted  by  the 
source  spectrum. 

From  Parse val’s  theorem,  the  mean  scattered  energy 

21 

over  the  resolution  footprint  is  proportional  to" 

E(rJr,r0)  =  J  \Q(f)\2(\<&s(rs\r,r0,f)\2)df,  (13) 

where  (|<l>J(rJ| r,r0,/)|2)  is  obtained  from  Eq.  (10). 

Sea  bottom  reverberation  is  typically  dominated  by  the 
incoherent  intensity.  In  this  case,  only  the  variance  term  in 
Eq.  (11)  contributes  to  the  total  intensity  in  the  right-hand 
side  of  Eq.  (13). 

III.  APPLICATION  TO  SEAFLOOR  REVERBERATION 
WITH  CALIBRATION 

In  this  section,  we  apply  the  scattering  model  developed 
in  Sec.  II  to  seafloor  reverberation.  Geophysical  parameters 
of  the  model  are  calibrated  with  measured  data. 

A.  OAWRS  images  of  the  ocean  environment 

Here  we  use  OAWRS  images  obtained  during  the  ONR- 
sponsored  Main  Acoustic  Experiment  20031  on  the  U.S.  East 
Coast  continental  shelf.  A  bistatic  OAWRS  system  compris¬ 
ing  of  a  vertical  source  array  and  a  horizontally  towed  re¬ 
ceiving  array  was  deployed  in  the  New  Jersey  Strataform 
located  200  km  south  of  Long  Island.  The  source  array  ex¬ 
cited  acoustic  waveguides  modes  that  propagated  tens  of  ki¬ 
lometers  in  the  continental  shelf  environment.  Transmissions 
consisted  of  1  s  duration  Tukey-shaded  linear  frequency 
modulated  pulses  in  three  different  frequency  bands;  390- 
440,  875-975,  and  1250-1400  Hz.  The  transmissions  were 
repeated  every  50  s  over  the  approximately  83  min  duration 
it  took  to  tow  the  receiving  array  along  a  10  km  track  line  at 
constant  tow  speed  of  2  m/s.  Data  were  collected  over  5-6 
tracks  each  day.  Here  we  focus  our  analysis  on  data  in  the 
390-440  Hz  band  as  this  was  the  wave  form  that  was  most 
frequently  used  for  imaging. 

Instantaneous  OAWRS  images  of  the  ocean  environ¬ 
ment  over  wide  areas  spanning  thousands  of  square  kilome¬ 
ters  are  formed  from  each  transmission.  The  arrivals  are 
charted  in  range  by  two-way  travel  time  and  matched  filter 
analysis,  and  in  bearing  by  beamforming.  ’  Scattered  fields 
from  objects  at  identical  horizontal  position  but  different 
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(a)  OAWRS  Data,  Normalized  Pressure  Level 


(b)  Simulated  Reverberation 


/  i 

'  / 

/'  / 

80  mV  /  , 

I  / 

/  # 

.  100  m'  f 

/  : 

/  ! 

\  r  - 

»  i  7 

'  120  m 


-30  -20 


-10  0  10 
Distance  from  source  (km) 


n 

tS  -10 
5 


-20 


-30 

-30 


-20  -10  0  10 
Distance  from  source  (km) 


FIG.  1.  (Color  online)  (a)  Instantaneous  OAWRS  image  of  the  ocean  environment  on  the  New  Jersey  continental  shelf.  The  image  is  averaged  over  five 
consecutive  transmissions  centered  at  09:32  EDT,  14  May  2003.  The  imaging  diameter  of  the  system  was  set  to  roughly  60  km  (40  s  imaging  time).  The 
moored  source  (the  white  diamond)  operating  at  390-440  Hz  is  the  coordinate  origin  at  39.0563N,  73.0365W.  The  towed  horizontal  receiving  array  moving 
along  51°E  heading  (indicated  by  black  line)  has  2.6°  azimuthal  resolution  at  array  broadside  and  15  m  range  resolution.  The  positive  vertical  axis  points 
north.  Depth  contours  are  indicated  by  dashed  lines.  The  image  has  been  corrected  for  source  level,  (b)  The  simulated  bottom  reverberation  corresponding  to 
the  data  in  (a),  calculated  using  the  computationally  efficient  approach  at  the  center  frequency  415  Hz  with  a  typical  observed  sound  speed  profile.  The  small 
black  1  km  by  1  km  box  shown  near  the  upper  left-hand  corner  indicates  the  location  used  for  calibrating  the  reverberation  model.  This  is  the  area  shown  in 
Fig.  6. 


depths  are  integrated  by  the  system  in  the  resulting  two- 
dimensional  (2D)  horizontal  OAWRS  images.  The  inherent 
left-right  ambiguity  in  charting  returns  with  a  horizontal  line 
array  receiver  is  evident  in  some  of  the  OAWRS  images  that 
will  be  presented  here. 

The  resolution  footprint  of  the  OAWRS  system  is  a 
function  of  both  range  and  azimuth.  The  receiving  array  has 
a  cross-range  resolution  given  by  p/3(<p)  where  the  angular 
resolution  f3(cf>)  is  given  by  /3((f>)  ~  1.44(\/L  cos  <^>)  for 
broadside  ((/>=())  through  angles  near  endfire  (<f>=  tt/2), 
where  A.  is  the  wavelength  of  the  signal  and  L  is  the  length  of 
the  array.  At  endfire,  the  angular  resolution  is  given  by 
/8(0=7r/2)~2.8yA./L.  A  Hanning  spatial  window  was  ap¬ 
plied  to  the  receiving  array  in  the  beamforming  to  reduce  the 

levels  of  the  sidelobes.  The  range  resolution  is  Ap=cl2B 

12 

after  matched  filtering,  where  B  is  the  signal  bandwidth. 

The  OAWRS  images  show  scattered  returns  from  large 
and  densely  populated  fish  groups,  small  or  diffuse  groups  of 
fish,  some  discrete  geologic  features  of  the  seafloor,  as  well 
as  continuous  diffuse  returns  from  the  seafloor.  A  typical 
image  is  shown  in  Fig.  1(a).1  Trends  in  the  images  include 
the  decay  of  scattered  intensity  with  increasing  bistatic  range 
due  to  acoustic  spreading  loss  and  attenuation  in  the  ocean 
waveguide,  and  spatial  variation  in  cross-range  resolution 
which  leads  to  spatially  varying  blurring. 
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B.  Estimating  scattering  level  of  the  environment 
from  OAWRS  data 

A  challenge  in  analyzing  sonar  imagery  is  to  distinguish 
diffuse  background  reverberation  from  discrete  returns  due 
to  targets  or  clutter.  The  discrete  returns  may  arise  from  sur¬ 
face  or  volume  inhomogeneities  such  as  marine  creatures, 
geologic  features,  or  man-made  objects.  It  is  often  not  known 
in  advance  from  where  in  the  water  column  the  targets  or 
clutter  features  arise. 

To  distinguish  diffuse  background  reverberation  from 
clutter  and  to  compare  the  scattering  level  of  returns  from 
different  ranges,  we  devise  a  two  step  approach.  First,  we 
detrend  an  OAWRS  intensity  image  by  correcting  for  mean 
transmission  over  the  water-column,  source  level,  and  the 
spatially  dependent  resolution  footprint  of  the  OAWRS  sys¬ 
tem.  This  leads  to  an  estimate  of  what  we  call  the  environ¬ 
mental  scattering  level  (ESL)  at  a  given  horizontal  location. 
Next  we  compute  statistics  of  ESL  variations  over  space  and 
time,  and  use  these  to  discriminate  diffuse  background  from 
clutter.  We  then  investigate  either  the  diffuse  background  or 
the  clutter.  Here  we  focus  on  the  background. 

Let  Lp(p^  |  r ,  r0 ,  t)  be  101og10  of  the  OAWRS  pressure 
squared  matched  filtered  output  centered  at  horizontal  loca¬ 
tion  ps  and  time  r,  which  may  include  some  spatial  and  tem¬ 
poral  averaging  in  pressure  squared  to  reduce  the  variance. 
Then  the  ESL(ps ,  r)  at  the  given  location  and  time  is  obtained 
by 
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Distance  from  source  (km)  Distance  from  source  (km) 

FIG.  2.  (Color  online)  (a)  A  zoomed-in  version  of  the  OAWRS  image  shown  in  Fig.  1(a)  in  the  south-east  region  containing  a  massive  fish  shoal,  (b) 
Corresponding  ESL  image.  Higher  ESL  levels  are  obtained  in  the  region  of  dense  fish  shoal  and  lower  levels  correspond  to  the  background.  The  transmission 
loss  used  for  correcting  this  image  is  shown  in  Fig.  3(c). 


ESL(p„?)  =  Lp(ps|r,r0,f)  -  SL(r0)  -  TL(p.s|r0) 

-TL(r|p,)-101ogIOA(p,),  (14) 

where  TL(ps|r0)  and  TL(r|ps)  are  101og10  of  the  average  of 
the  antilog  of  transmission  loss  over  the  water  column.  The 
former  is  for  transmission  from  source  location  to  environ¬ 
mental  location  ps  and  the  latter  is  from  this  environmental 
location  to  the  receiver.  The  spatially  varying  resolution  foot¬ 
print  of  the  receiving  array  has  area  A(ps)  =  [/3((f>)  /  2tt][tt(p 
+  Ap)2- 7rp2]  —  p/3{<f>)Ap,  and  SL(r0)  is  the  source  level.  En¬ 
vironmental  scattering  level  is  then  10  log10  of  a  ratio  and 
has  units  simply  of  deciBels,  like  scattering  strength. 

Here  we  use  a  statistical  propagation  model  based  on  the 
range-dependent  parabolic  equation  calibrated""  against 
one-way  transmission  data  from  the  same  set  of  measure¬ 
ments  to  estimate  the  mean  transmission  loss  over  the  water- 
column  and  obtain  TL(ps|r0)  and  TL(r|ptj.  Effects  such  as 
internal  waves,  eddies  and  turbulence  cause  fluctuations  in 
ocean  sound  speed  and  density  structure.  These  lead  to  scin¬ 
tillation  in  ocean  acoustic  waveguide  propagation  caused  by 
changes  in  the  multimodal  interference  structure.  In  this  en¬ 
vironment,  the  acoustic  waveguide  modes  were  found  ex¬ 
perimentally  to  become  sufficiently  randomized  beyond 
roughly  2  km  in  range  from  the  source""  that  the  expected 
acoustic  intensity  varies  only  gradually  with  increasing  range 
and  does  not  exhibit  the  strong  periodic  modal  interference 
structure  in  range  and  depth  present  in  nonrandom 
waveguides.16'"0  The  antilog  of  the  mean  transmission, 
TL(pj|r0)  and  TL(r|pJ  are  proportional  to  the  expected 
acoustic  power  transmitted  through  the  waveguide  at  any 
given  range  and  azimuth. 

Figures  2(a)  and  2(b)  are  OAWRS  images  of  the  mean 
matched  filtered  scattered  intensity  and  the  corresponding 
ESL  of  the  ocean  environment  on  14  May  2003  after  inco¬ 
herently  averaging  intensity  over  five  consecutive  transmis¬ 
sions  (pings)  centered  at  time  instance  9:32  Eastern  Daylight 
Time  (EDT).  The  ESL  image  in  Fig.  2(b)  is  derived  from  Eq. 
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(14)  where  images  corresponding  to  the  right-hand  side  of 
Eq.  (14)  appear  in  Figs.  2(a)  and  3.  The  mean  transmission 
losses  shown  in  Figs.  3(a)  and  3(b)  are  in  good  agreement 
with  measurements  or  models  of  transmission  losses  reported 
in  Refs.  22-24  for  the  New  Jersey  Strataform. 

In  order  to  distinguish  the  statistically  stationary  back¬ 
ground  ESL  from  that  due  to  moving  biological  clutter,  we 
average  N  consecutive  OAWRS  images,  which  already  were 
formed  from  five-ping  running  averages,  to  obtain  a  mean 
OAWRS  image  and  a  standard  deviation  OAWRS  image 
over  a  wide  spatial  area.  Figures  4(a)  and  4(b)  show  the 
mean  and  standard  deviation  of  Af=361  ESL  images  derived 
from  data  acquired  within  the  frequency  band  from 
390  to  440  Hz  on  14  May  2003.25 

We  find  in  our  2003  OAWRS  images  that  (a)  high  mean 
and  high  standard  deviation  characterizes  regions  where  fish 
are  regularly  passing  through,  (b)  high  mean  and  low  stan¬ 
dard  deviation  characterizes  regions  where  fish  consistently 
remain,  (c)  low  mean  and  high  standard  deviation  character¬ 
izes  regions  where  fish  pass  through  intermittently,  and  (d) 
low  mean  and  low  standard  deviation  characterizes  regions 
of  diffuse  reverberation  from  the  seafloor.  This  can  be  seen  in 
Fig.  5(a)  where  ESL  is  plotted  as  a  function  of  sequential 
ping  number  for  various  spatial  locations  from  Fig.  4(b).  In 
particular,  ESL  in  regions  of  dense  fish  population  is  roughly 
20-25  dB  higher  than  in  regions  where  diffuse  seafloor  re¬ 
verberation  is  the  dominant  scattering  mechanism.  Note  that 
the  above-mentioned  (b)  could  also  be  caused  by  static  geo¬ 
logic  features  but  was  ruled  out  in  this  experiment  as  data 
from  the  area  shown  in  Fig.  4,  six  days  prior,  showed  the 
scattered  levels  returning  to  those  of  the  background  when 
the  massive  fish  shoal  was  absent. 

For  estimating  the  background  ESL  due  to  the  seafloor, 
we  focus  analysis  on  regions  with  upslope  bathymetry  rela¬ 
tive  to  the  incident  wave  that  had  ambiguous  regions  that 
were  downslope.  This  is  to  help  eliminate  the  effect  of  left- 
right  ambiguity,  as  seafloor  scattering  has  been  shown  to  be 
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FIG.  3.  (a)  TL(pt|r0)  and  (b)  TL(r|p,)  are  10  log10  of  the  mean  antilog  transmission  losses  over  the  water  column  from  source  and  receiver  respectively  to 
the  area  shown  for  the  390-440  Hz  band,  (c)  The  quantities  in  (a)  and  (b)  blurred  with  the  spatially  dependent  resolution  footprint  of  the  receiving  array. 


more  intense  when  bathymetry  is  upward  sloping  in  the  di¬ 
rection  of  the  incident  wave,  and  weaker  when  bathymetry  is 
downward  sloping,  13,14'~6~~ 9  as  expected  from  physical  prin¬ 
ciples.  Temporal  variations  in  ESL  levels  for  various  spatial 
locations  with  upslope  bathymetry  are  shown  in  Fig.  5(b). 
The  area  shown  in  Fig.  4  consists  of  down  sloping  bathym¬ 
etry  but  is  shown  to  enable  comparison  of  background  ESL 
levels  with  those  from  fish  shoal. 

The  mean  seafloor  ESL  obtained  here  (and  tabulated  in 
the  forthcoming  Table  II)  can  now  be  used  to  calibrate  the 
reverberation  model  presented  in  Sec.  II  for  seafloor  scatter¬ 
ing  in  the  Strataform  area.  The  standard  deviation  of  ESL  for 
seafloor  scattering  has  been  found  by  temporal  and  spatial 
averaging  to  be  about  2  dB.  This  standard  deviation  will  be 
useful  in  determining  the  accuracy  of  parameter  estimates 
needed  to  calibrate  the  reverberation  model. 

C.  Estimating  geologic  parameters  necessary  for 
reverberation  model 

The  geological  parameters  necessary  to  implement  the 
reverberation  model  defined  in  Eqs.  (7)  and  (11)  are  esti¬ 
mated  in  the  Appendix.  In  particular,  the  mean,  standard  de- 
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viation,  and  correlation  of  T*  and  T  d  are  estimated  from  the 
statistical  moments  of  sound  speed  and  density  provided  in 
Refs.  30,  31,  and  3.  The  results,  tabulated  in  Table  I,  are 


TABLE  I.  Statistical  geologic  properties  of  the  New  Jersey  Strataform  used 
to  calibrate  the  bottom  reverberation  model. 


Statistical  parameters 

Average  values  over  the  region 

c  (m/s) 

1700 

<jc  (m/s) 

35.5 

fj.2=crclc 2 

4.6889e-04 

d  (g/cm3) 

1.90 

<rd  (g/cm3) 

0.16 

rf=<?d / 

6.2e-3 

0L> 

0 

0.0083 

VardL) 

0.0083 

<r,> 

-0.0063 

<rs> 

0.0066 

V  ar(r  d) 

0.0065 

Cov(TK,Td) 

-0.0065 
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Distance  from  source  (km)  Distance  from  source  (km) 

FIG.  4.  Temporal  and  spatial  variations  in  ESL  for  the  area  in  Fig.  2  for  data  in  the  390-440  Hz  band  collected  on  14  May  2003.  (a)  Mean  and  (b)  standard 
deviation  of  361  ESL  images  formed  from  transmissions  over  a  measurement  time  period  ranging  from  09:10  to  17:20  EDT. 


consistent  with  those  found  in  Ref.  7  which  are  based  on 
measurements  from  ten  different  shallow  and  deep  water 
sites  around  the  world. 

From  Table  I  we  observe  that  (I^)2  and  (rd)2  are  at  least 
2  orders  of  magnitude  smaller  than  Var(T J  and  Var(rd). 
The  coherent  term  in  Eq.  (10)  will  therefore  be  negligible  in 
comparison  to  the  incoherent  term.  The  total  reverberation 
intensity  can  then  be  approximated  as  the  variance  of  the 
scattered  held  using  Eqs.  (13)  and  (11). 

D.  Maximum  likelihood  estimator  for  coherence 
volume  Vc 

To  implement  the  reveberation  model  in  Eqs.  (13)  and 
(11)  at  a  given  location,  an  estimate  is  needed  for  the  coher¬ 
ence  volume  Vc  of  the  random  seabed  inhomogeneities.  An 


expression  for  the  maximum  likelihood  estimator  (MLE)  for 
the  coherence  volume  Vc  in  terms  of  the  measured  reverbera¬ 
tion  data  and  model  is  provided  here. 

The  instantaneous  scattered  held  measured  at  the  re¬ 
ceiver  can  be  modeled  as  a  circular  complex  Gaussian  ran- 
dom  variable.  The  log  transform  of  the  matched  hltered 
intensity  within  the  resolution  footprint  centered  at  ps  then 
obeys  a  Gaussian  probability  density  when  the  time- 
bandwidth  product  is  sufficiently  large,  as  does  the  ESL, 

p(ESUpM) 

=  (^)  exp(-  ^[ESL(p.s)  -<ESL(Pj,Vf)>]2), 

(15) 

where  [jl  is  the  time-bandwidth  product,  which  for  the  cur- 


FIG.  5.  (a)  Temporal  variation  in  ESL  at  various  spatial  locations  as  a  function  of  transmission  number  or  ping,  for  the  image  in  Fig.  4.  Gray  lines  in  (a)  are 
ESL  at  locations  with  persistently  strong  scattered  returns  from  fish  shoals.  Black  lines  indicate  ESL  at  locations  consistent  with  the  background,  (b)  Temporal 
variations  in  ESL  for  the  region  with  upslope  bathymetry  in  the  vicinity  of  the  small  box  shown  in  Fig.  1(b). 
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rent  processing  is  greater  than  10.  This  is  because  the 
matched  filtered  intensity  data  are  averaged  over  two  inde¬ 
pendent  resolution  cells  and  five  independent  pings. 

As  /x  is  large,  the  bias  due  to  log  transformation,  defined 
in  Eq.  (22)  of  Ref.  32,  is  negligible.  We  can  then  make  the 
approximation 

<ESL(ps,Vc)>  =  10  log10(W(ps|r,r0))  +  10  log10(Vc) 

-<TL(pJ|r0))-<TL(r|pi)),  (16) 

where  we  have  introduced  the  new  variable 
10  logl0(W(pv|r,r0)) 

=  10  log10[£(pjr,r0)]  -  10  log10(Vc),  (17) 

because  it  does  not  depend  on  Vc.  Here,  E( pv  |  r ,  r0)  is  defined 
in  Eq.  (13),  and  is  proportional  to  the  total  energy  returned 
from  the  resolution  footprint  of  range  extent  A r.  It  is  advan¬ 
tageous  to  express  our  intensity  data  in  logarithmic  units  as  it 

transforms  signal-dependent  noise  into  additive  signal- 

32 

independent  noise. 

The  ESL  of  the  bottom  measured  in  N  different  reso¬ 
lution  footprints  of  OAWRS  imagery  for  s  =  1 , 2,3, ...  ,N  can 
be  assumed  to  be  independent  and  identically  distributed. 
Their  joint  probability  density  function  is 

p(ESL(p,),ESL(p2),  ...  ,ESL(pN)\Vc) 

=  II  (  ^)  exp(-  ^[ESL(p.s)  -  <ESL(ps,  Vc)>]2)  . 

(18) 

The  maximum  likelihood  estimator  Vc  for  the  coherence  vol¬ 
ume  Vc  is  given  by 

-^r{ln[p(ESL(p1),ESL(p2),  . . .  ,ESL(pw)|Vc)]} 
dVc  vc=Vc 

=  0.  (19) 

Expanding  Eq.  (19)  leads  to 

ioiog10(yc) 

1  v  ,  -v 

=  “2  ESL(p.5)  -  -2  {10  log|0(W/(r|r0,ps)) 

Jy  s=l  Jy  S=1 

-TL(pJr0)»TL(r|ps)}  (20) 

the  desired  estimator. 

E.  Reverberation  model  calibration 

Implementation  of  the  full  reverberation  model  in  Eq. 
(11)  for  range-dependent  waveguides  involves  numerically 
calculating  3D  spatial  derivatives  of  the  complex  waveguide 
Green  function.  This  leads  to  large  3D  matrices  that  makes 


FIG.  6.  Full  bistatic,  390-440  Hz  broadband,  depth-integrated  bottom  re¬ 
verberation  model  simulation  for  a  horizontal  transect  across  the  1  km  by 
1  km  area  shown  in  Fig.  1(b).  The  monopole  term,  dipole  term,  and  cross¬ 
term  contributions  to  the  reverberation  intensity  are  plotted  for  comparison. 


the  implementation  cumbersome  and  reduces  computational 
speed.  For  instance  at  400  Hz,  the  spatial  sampling  required 
is  approximately  0.2  m  in  range,  0.5  m  in  depth,  and  roughly 
3°  in  azimuth  for  computational  accuracy. 

To  avoid  this  problem,  we  calibrate  the  full  reverbera¬ 
tion  model  using  only  a  small  area,  1  km  by  1  km,  shown  in 
Fig.  1(b).  The  coherence  volume  Vc  is  then  determined  with 
the  approach  outlined  in  Sec.  Ill  D.  We  provide  a  method  for 
implementing  the  reverberation  model  that  is  computation¬ 
ally  efficient  over  wide  areas  in  the  next  section. 

The  scattered  intensity  of  Eq.  (11)  is  a  sum  of  (1)  a 
purely  monopole  term  dependent  on  the  variance  of  frac¬ 
tional  change  in  compressibility,  (2)  a  purely  dipole  term 
dependent  on  the  variance  of  fractional  change  in  density, 
and  (3)  the  covariance  between  the  two.  The  contribution  of 
these  three  terms  for  a  horizontal  transect  through  the  small 
area  is  shown  in  Fig.  6. 

The  MLE  of  the  coherence  volume  Vc  is  summarized  in 
Table  II.  We  assume  the  inhomogeneites  are  isotropic,  so  that 
Vc  =  ^tt12c,  where  lc  is  the  coherence  radius.  The  coherence 
lengths  in  x,  y,  and  z  directions  are  lx=ly=lz=2lc.  The  acous¬ 
tically  determined  correlation  length  obtained  here  is  compa¬ 
rable  to  the  geologic  correlation  depth  of  between  0.3  and 

33 

3  m  from  analysis  of  core  data. 

As  the  background  ESL  has  a  standard  deviation  of 
±2  dB,  the  estimate  101og10  V r  also  has  the  same  standard 
deviation.  This  leads  to  a  correlation  length  in  the  range  from 
15  to  21  cm.  There  could  also  be  uncertainties  in  our  esti¬ 
mates  of  the  moments  of  F(<.  and  T d.  The  variations  in 
Var(T  J  needed  to  match  the  data  over  the  range  of  uncer¬ 
tainty  in  the  correlation  length  occupies  roughly  a  factor  of  2. 


TABLE  II.  Estimates  for  the  New  Jersey  Strataform.  The  standard  deviation  for  ESL  and  SS  is  ±2  dB. 


/(Hz) 

ESL  (dB) 

Vc  (m3) 

2lc  (m) 

fd  (dB) 

fP  (m) 

SS  (dB) 

390-440 

-67.5 

0.0030 

0.18 

4.2 

0.46 

-37 
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Monostatic  reverberation  in  Pekeris  waveguides  (390-440  Hz) 


Range  (km) 


Range  (km) 


FIG.  7.  Broadband  390-440  Hz,  monostatic  reverberation  in  Pekeris  waveguides  with  (a)  sand  and  (b)  silt  bottom  for  a  0  dB  re  1  /zPa  at  1  m  source  level 
comprising  of  a  single  element  at  50  m  depth.  The  receiving  array  parameters  are  the  same  as  those  used  in  the  experiment  in  Fig.  1(a).  The  Green  function 
for  Pekeris  waveguides  calculated  using  the  normal  mode  code,  Kraken  (Ref.  36).  (c)  The  dipole  scaling  factor  equal  to  the  ratio  between  the  full  term  and 
the  monopole  term  plotted  in  deciBels  remain  fairly  constant  to  within  ±0.1  dB  and  is  independent  of  the  bottom  type,  (d)  Comparison  of  the  magnitudes  of 
the  Greens  function  and  its  spatial  derivative  at  the  water-sediment  interface  at  415  Hz. 


Here  we  ignored  a  potential  degradation  due  to  matched  fil¬ 
tering  the  data  in  an  ocean  waveguide  while  calibrating  them 
with  the  reverberation  model.  This  degradation  is  typically 

22  34 

small,  less  than  3  dB  for  the  ranges  of  interest.  ’ 

IV.  COMPUTATIONALLY  EFFICIENT  IMPLEMENTATION 

When  working  on  a  field  experiment,  it  is  often  instruc¬ 
tive  to  rapidly  predict  seabed  reverberation  over  wide  areas, 
to  aid  in  experimental  design.  A  computationally  efficient 
approach  for  computing  seabed  reverberation  in  range  de¬ 
pendent  waveguides  over  wide  areas  is  developed  here.  It  is 
based  on  the  introduction  of  a  dipole  scaling  factor  and  a 
depth-penetration  factor  that  simplify  the  numerical  imple¬ 
mentation. 

We  find  it  a  good  approximation  to  take  the  dipole  and 
cross  terms  in  the  reverberant  intensity  of  Eq.  (11)  to  be 
proportional  to  the  monopole  term.  This  is  illustrated  for 
horizontally  stratified  Pekeris  sand  and  silt  waveguides  in 
Figs.  7(a)  and  7(b)  with  collocated  source  and  receiver,  and 
for  the  bistatic  scenario  in  a  range  dependent  waveguide  in 
Fig.  6.  The  total  intensity  can  be  estimated  from  the  mono- 
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pole  term  alone  by  multiplication  with  a  constant  factor,  Fd, 
which  we  call  the  dipole  scaling  factor.  This  can  be  seen  by 
noting  in  Fig.  7(c)  that  the  ratio  of  the  total  intensity  to  the 
monopole  term  is  a  constant,  with  10  log10T’d~4.2  dB.  This 
approximation  is  a  consequence  of  the  fact  that  |VG|  is  pro¬ 
portional  to  k\G\,  as  shown  in  Fig.  7(c),  which  means  that  the 
dipole  and  monopole  terms  are  proportional.  Without  this 
approximation,  it  is  currently  impractical  to  compute  rever¬ 
beration  over  wide  areas  because  of  the  tremendous  compu¬ 
tational  effort  needed  to  determine  the  derivatives  of  the 
Green  function  in  the  dipole  term. 

The  depth  integral  of  Eq.  (11)  for  volume  scattering  can 
be  approximated  by  the  seafloor  surface  contribution  multi¬ 
plied  by  a  depth  penetration  factor,  Fp.  The  depth  penetration 
factor  provides  the  equivalent  depth  over  which  acoustic  in¬ 
tensity  in  the  seabed  is  of  similar  magnitude  as  that  on  the 
seafloor,  before  it  decays  rapidly  with  depth  due  to  sediment 
attenuation.  This  factor  is  simply  the  ratio  of  the  volume 
integral  to  the  seafloor  contribution,  which  are  shown  in 
Figs.  8(a)  and  8(b),  with  corresponding  Fp  ratio  in  Fig.  8(c), 
which  is  found  to  be  dependent  on  bottom  type,  roughly 
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Monostatic  reverberation  in  Pekeris  waveguide  (390-440  Hz) 


FIG.  8.  Comparison  of  the  390-440  Hz  broadband  depth-integrated  bottom  reverberation  intensity  with  the  broadband  intensity  at  the  water-sediment 
interface  in  the  Pekeris  (a)  sand  and  (b)  silt  waveguides,  (c)  The  effective  penetration  depth  in  meters  plotted  for  Pekeris  sand  and  silt  waveguides. 


0.5  m  for  sand  and  1  m  for  silt.  The  depth  penetration  factor 
is  also  a  function  of  frequency,  decreasing  with  increasing 
frequency. 

Another  approximation  that  can  greatly  increase  compu¬ 
tational  efficiency  is  to  calculate  the  scattered  field  at  just  the 
center  frequency  of  the  broadband  signal.  Figure  9  compares 
the  full  broadband  depth  integrated  monostatic  reveberation 
to  the  range-averaged  scattered  field  at  the  center  frequency 
with  the  two  approximations  introduced  earlier  applied  to  the 
result.  The  latter  single  frequency  range-averaged  approxi¬ 
mation  provides  a  good  estimate  to  the  broadband  result. 

The  approximations  made  in  this  section  require  only 
the  Green  function  magnitude  on  the  seabed  surface  which 
can  be  stored  in  sparse  2D  matrices  rather  than  dense  3D 
matrices  required  for  the  full  calculation.  The  reverberation 
intensity  is  accurate  to  within  ±0.5  dB  for  ranges  greater 
than  3  km  and  to  within  3  dB  for  ranges  between  1  and 
3  km.  Figure  1(b)  illustrates  this  approach  applied  to  model 
reverberation  over  wide  areas  on  the  New  Jersey  continental 
shelf  corresponding  to  the  OAWRS  data  Fig.  1(a).  The  rever¬ 
beration  levels  from  the  simulation  are  comparable  to  the 
background  levels  in  the  data,  especially  in  areas  adjacent  to 
the  massive  fish  shoal. 
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V.  SCATTERING  STRENGTH 

Here  we  derive  the  scattering  strength  (SS)  of  the  bot¬ 
tom  on  the  New  Jersey  continental  shelf  using  the  model  and 
calibration  results.  As  the  seabed  inhomogeneities  have  co- 
herence  lengths  that  are  smaller  than  the  wavelength,"  we 
can  apply  the  sonar  equation  to  approximate  SS  of  the  bot¬ 
tom  from  Eq.  (11)  with  the  approximations  introduced  in 
Sec.  IV  as 

SS  ~  10  logl0[FrfFpF(±4Var(r,)] .  (21) 

The  SS  of  the  sandy  bottom  on  the  New  Jersey  Strataform  is 
tabulated  in  Table  II.  Analysis  with  a  normal  mode  propaga¬ 
tion  model3'1  indicates  that  the  first  10  modes  are  important 
contributors  to  the  acoustic  field  at  a  range  of  20  km  from 
the  source  where  the  calibration  was  done.  At  this  range,  the 
depth-averaged  intensity  of  mode  10  is  about  1/5  the  depth- 
averaged  intensity  of  mode  1 .  The  grazing  angles  for  modes 
1-10  range  from  2°  to  12°  and  the  derived  SS  level  of 
roughly  -37  ±  2  dB  are  determined  by  these  grazing  angles 
after  long  range  propagation  in  the  waveguide  with  the  effect 
of  multimodal  coupling  in  the  scattering  process. 
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FIG.  9.  Full  390-440  Hz  broadband  depth-integrated  monostatic  bottom 
reverberation  compared  to  simulated  result  at  center  frequency  415  Hz  cal¬ 
culated  using  efficient  approach  with  the  dipole  scaling  and  depth  penetra¬ 
tion  factors. 


For  reverberation  measured  at  390-440  Hz  bandwidth 
considered  in  this  paper,  volume  scattering  is  expected  to  be 
the  dominant  mechanism  at  all  grazing  angles  for  a  sandy 
seabed  as  shown  in  Figs.  10  and  11  of  Ref.  5.  Our  estimate 
agrees  well  with  those  of  other  ocean  environments  with 
sandy  bottoms  in  Refs.  5  and  10,  suggesting  low  surface 
roughness  at  this  site  and  the  consequent  importance  of  eva¬ 
nescent  energy  in  generating  low-angle  volume  backscatter 
via  subcritical  ensonification. 

The  overall  frequency  dependence  of  SS  in  Eq.  (21)  is 
determined  by  the  wave  number  k,  which  is  directly  propor¬ 
tional  to  frequency,  the  depth  penetration  factor  fp  that  de¬ 
creases  with  increasing  frequency,  and  potentially  the  coher¬ 
ence  volume  sampled  by  the  acoustic  wave  Vc.  The 
coherence  volume  may  depend  on  frequency  as  short  wave¬ 
length  acoustic  fields  are  more  sensitive  to  small  scale  fluc¬ 
tuations  in  medium  properties  than  long  wavelength  acoustic 
waves.  This  acoustically  determined  coherence  volume  may 
differ  from  the  geologic  coherence  volume  that  is  determined 
purely  from  spatial  variations  in  medium  properties. 

VI.  CONCLUSION 

An  analytic  model  for  3D,  bistatic  scattering  from  me¬ 
dium  inhomogeneities  is  developed  from  first  principles  by 
application  of  Green’s  theorem.  Statistical  moments  of  the 
scattered  field  are  expressed  in  terms  of  statistical  moments 
of  medium  compressibility  and  density  fluctuations.  The 
model  is  applied  to  seabed  reverberation  and  optimally  cali¬ 
brated  with  both  OAWRS  and  geological  data  on  the  New 
Jersey  continental  shelf.  Analysis  with  the  model  indicates 
that  (1)  seabed  reverberation  is  incoherent,  and  (2)  scattering 
strength  varies  with  frequency  depending  on  wavenumber  k, 
medium  coherence  volume  Vc  and  seabed  depth  penetration 
factor  Fp  following  a  10  log l0(FpVck4)  dependence. 

An  efficient  numerical  approach  is  also  developed  for 
rapidly  computing  seabed  reverberation  over  wide  areas  for 
bistatic  sonar  systems  in  range-dependent  ocean  waveguides. 


It  exploits  the  correlation  between  monopole  and  dipole  scat¬ 
tering  terms  and  the  limited  penetration  of  acoustic  fields  in 
the  seabed.  The  model  handles  the  scattering  of  evanescent 
waves  in  the  seabed.  Finally,  an  approach  for  distinguishing 
the  statistically  stationary  background  reverberation  from  the 
scattered  fields  of  moving  targets  in  sonar  data  by  tracking 
the  temporal  and  spatial  evolution  of  the  returns  is  presented. 
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APPENDIX:  DERIVATION  OF  STATISTICS  OF 
GEOACOUSTIC  PARAMETERS 


The  seabed  on  the  New  Jersey  Strataform  is  predomi¬ 
nantly  comprised  of  sand.  It  has  been  mapped  by  geophysi¬ 
cal  surveys  with  analysis  provided  in  Refs.  3,  30,  and  31. 
From  these,  we  gathered  the  statistical  moments  of  seabed 
sound  speed  and  density  fluctuations  and  list  them  in  Table  I. 
We  assume  that  sound  speed  N(c,crc)  and  density 

d,(pr)  ~  N(d,  <id)  obey  Gaussian  random  processes  that  are 
locally  stationary  in  a  given  area.  These  statistical  moments 
can  vary  between  areas  to  account  for  potential  changes  in 
sediment  type.  We  assume  that  sound  speed  and  density  are 
uncorrelated,  which  is  a  valid  assumption  given  the  large 
scatter  in  the  data  in  Fig.  10(b)  of  Ref.  30.  Despite  this  as¬ 
sumption,  rK  and  r(/  can  still  be  correlated. 

From  Eqs.  (2)  and  (3),  the  means  are  given  by 

<r*(rf)>  =  0  (Al) 

and 

(A2) 


where  d  and  c  as  defined  before  are  the  mean  density  and 
sound  speed  within  the  sonar  resolution  footprint,  respec¬ 
tively.  The  second  moments  are  given  by 


<r3(r,)> 


(A3) 


and 


<r2K(rf)>: 


(A4) 


the  variances  by 

Var(rrf(r,))  =  (r2(r/))-(r,(r,))2  (A5) 


and 

Var(rK(rf))  =  <r2(r,)>  -  <r>,)>2,  (A6) 

and  their  covariance  by 
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cov(rK,rrf)  =  \  l- 


d,(  r,) 


K,(  r,) 


-  1 


d,{ r,) /  k\ d*(rt)c*(rt) /  ' 


(A7) 


where 


v  d,(rt)  /  \  c~{ r,) , 


(A8) 


Table  I  provides  a  summary  of  these  calculations  for  the 
New  Jersey  Strataform. 
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